Avoiding dissipation in a system of three quantum harmonic oscillators 
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We analyze the symmetries in an open quantum system composed by three coupled and detuned 
harmonic oscillators in the presence of a common heat bath. It is shown analytically how to engineer 
the couplings and frequencies of the system so as to have several degrees of freedom unaffected by 
decoherence, irrespective of the specific spectral density or initial state of the bath. This partial 
thermalization allows observing asymptotic entanglement at moderate temperatures, even in the 
nonresonant case. This latter feature cannot be seen in the simpler situation of only two oscillators, 
highlighting the richer structural variety of the three-body case. When departing from the strict 
conditions for partial thermalization, a hierarchical structure of dissipation rates for the normal 
modes is observed, leading to a long transient where quantum correlations such as the quantum 
discord are largely preserved, as well as to synchronous dynamics of the oscillators quadratures. 

PACS numbers: 



I. INTRODUCTION 



Prevention of decoherence and dissipation in open 
quantum systems is a fundamental condition for the pres- 
ence of quantum phenomena in warm macroscopic every- 
day world. While decoherence has been extensively stud- 
ied from the early 1980s Ju_2[ to nowadays (for a list of 
reviews see for instance |3|-l6|), and different mechanisms 
to avoid it have been discussed in the last years, including 
strategies to engineer it for applications 0-12I- Further- 
more, some macroscopic systems from marine algae [10| 
to metal carboxylates [11], can present quantum correla- 
tions at high temperatures. These recent experiments 
suggest that avoiding a complete quantum-to-classical 
transition induced by monitoring of the environment is 
not only artificially possible but can occur inherently 
in natural phenomena. The mechanisms that produce 
such survival or even construction of coherence at large 
time scales remain almost unclear, but different theoret- 
ical strategies have been proposed in order to predict it, 
mostly motivated in the context of quantum computation 
p^2l-[l6j . Indeed decoherence and the leak of information 
to the environment are the major obstacle in quantum 
processing of information and construction of quantum 
memories [l7j . 

In this context, one of the strategies to bypass decoher- 
ence is exploiting dynamical symmetries in the system- 
environment interaction. In order to generate unitary 
evolution in a certain subspace of the Hilbert space of an 
open system, a common dissipation (where several units 
equally couple to the same environment) has been firstly 
used in a two-qubit s^tem ^18h21] and later extended to 
multiple qubits |22l , |23| and continuous variable systems 
[24J29I, I31I , I33J . A general framework has been devel- 
oped with several contributions (see for example [15| and 
the references therein) agglutinating the main concepts of 
decoherence- free subspaces (DFSs), Noiseless subsystems 



(NSs) or more recently information-preserving structures 
(IPSs) [l6| and proposing general methods to obtain 
them for arbitrary Hamiltonians. In our work we extend 
previous studies in the context of continuous variables ex- 
ploring the vaster landscape offered when moving from 
two to three coupled harmonic oscillators in the search 
for NSs induced by non-trivial symmetries in the system. 

Previous works on dissipative harmonic oscillators re- 
ported that in presence of identical frequencies and cou- 
plings between oscillators the symmetry of the collective 
motion and their interaction with the environment can 
lead to the effective decoupling of some of the normal 
modes of the system from the bath [24J28J . The cases of 
different frequencies [29, 30] or couplings [31| is instead 
less studied and understood. The natural step of con- 
sidering three harmonic oscillators beyond the symmet- 
ric configuration of identical oscillators already provides 
much more phenomenological richness while at the same 
time it allows for analytic treatment and gives valuable 
intuition when pursuing a further extension to the case 
of N oscillators. 

We show in this work how to obtain NSs and bypass 
decoherence independently of the bath properties such as 
temperature or frequency distribution, considering differ- 
ent frequencies, couplings and boundary conditions, in 
the presence of a common bath. We also analyze how 
by using these NSs, quantum correlations like entangle- 
ment can persist (by two different mechanisms) in the 
asymptotic limit of the dynamical evolution, given in our 
case by a Markovian-type dynamics in the weak-coupling 
limit. Furthermore even for three different oscillators an 
homogeneous dynamics, such as collective quantum syn- 
chronization, is shown to be present in some parameters 
manifolds. Synchronization is a widespread phenomenon 
in physical, biological, chemical and social systems p2l , 
only recently explored in the quantum regime [33|, [SJ . 
Our results for the system of three oscillators in presence 
of common dissipation may be implemented with ions 



in linear Paul traps by following the proposal in [35|. 
Experimental realization of coupled harmonic oscillators 
appear in optical [36|, [33 and superconducting [38|, [39| 
cavities as well as trapped ions [40|, |41| or nanoelectrome- 
chanical resonators \^. Three coupled elements archi- 
tectures are also known to allow for isochronous synchro- 
nization of semiconductor lasers with delayed coupling or 
neuronal models [43]. 

The paper is organized as follows: in Sect. [U we in- 
troduce the model for the system of three harmonic os- 
cillators dissipating into a common bath, also in terms of 
the normal modes of the system. For certain particular 
values of the system's parameters (and independently of 
the bath characteristics), one of several normal modes 
can be protected from decoherence. We find analytically 
these conditions leading for NSs in Sect. IIIII giving some 
specific cases that are analyzed in more detail. In the 
following we focus on the long time dynamics of the sys- 
tem considering in Sect. IIII Al the creation and conserva- 
tion of entanglement at asymptotic times. Section IIII Bl 
is dedicated to ability of the system to generate asymp- 
totic synchronized states and exploring its connection to 
the robustness of quantum correlations. In section HVl 
deviations of the NSs conditions are considered. This 
leads to dynamical relaxation of the system that con- 
verges towards a thermal state. We conclude with Sect. 
IVl summarizing our main results. 



II. THE MODEL 

We start with a Hamiltonian describing three coupled 
quantum harmonic oscillators with arbitrary frequencies 
and coupling constants between them. For simplicity we 
suppose unit masses and natural units {h = ks = I)'- 
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where Pi and Qi represents respectively the momen- 
tum and position operators of each harmonic oscillator 
{[Qi^Pj] = '^^ij) and A^ = 3 for our case. This equation is 
conveniently expressed in quadratic matrix form as: 
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where I is the identity {N x N) matrix, q^ = (g'l, ..., ^at) 
and H contains all the parameters of the system, i.e. 
the squared frequencies and couplings between oscilla- 
tors. We will only consider H with positive eigenvalues, 
so as to have bounded states (attractive potential). 

The environment is introduced by equally coupling 
each oscillator of the system to the same thermal bath, 
which is described by an infinite collection of independent 
bosonic modes: 
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where [Xc,n/3] = i6af3 and masses are unitaries. We 
will use along the paper greek subindices to refer to bath 
modes while latin ones are reserved for system oscillators 
(i, j) and normal modes (/c, n). The system-bath interac- 
tion reads: 



N 



Hi = 2_^ Qi /_^ A^Xq,, 



(4) 
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with a factorized form Ri = S®B of an operator S acting 
only on the system's Hilbert space and B acting on the 
environment one. As usual, this type of interaction yields 
a renormalization of the frequencies that we may include 

directly in our model by changing a;| -^ uof + ^^ 2#" 

The normal modes basis of the system ([2|) is obtained 
after a canonical transformation of the system operators 
through the orthogonal basis-change matrix J^\ 
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which diagonalizes 1-L {ci^l-Lci = Q^l^Q). Here Q = 
F^l-LF is a diagonal matrix containing the squared nor- 
mal modes frequencies Q.ri with n = 1,2,..., A". In this 
basis Hs now represents the Hamiltonian for a A' = 3 
uncoupled harmonic oscillators, or normal modes, related 
with the original (natural) modes by J^, so we can write: 



1 ^ 
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and the interaction Hamiltonian of Eq. (|4]) is now: 
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We stress that even if the oscillators are coupled with the 
same strength to the bath, the couplings of the normal 
modes to the bath, a^^, are not homogeneous: 
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The effective couplings tZn only depend on the canoni- 
cal transformation, i.e. on the system's parameters and 
arrangement defined by H. This suggests a strategy to 
protect one or more normal modes from the environment 
action based on proper tuning of {cj^, Xij}. Our analysis 
in section IIIII addresses this point while deviations form 
the condition of vanishing effective coupling of a system 
normal mode and the environment are explored in section 

We mention that while this work focuses on the case 
of three coupled harmonic oscillators, the description in 
terms of effective couplings here provided is rather gen- 
eral and applies for arbitrary networks of A' harmonic 



oscillators. The case of a common bath for all oscillators 
in the system corresponds to situations where the correla- 
tion length in the environment is larger than the system 
size. This assumption is not crucial for the framework 
we develop here, though any other choice would produce 
different specific analytic expressions. Finally, the case 
of a separate bath for each oscillator was shown in Ref. 
[33| to yield equal decoherence for all normal modes and 
therefore no NSs nor synchronization. 



A. Unbalanced bath couplings 

While the common bath situation is mostly a charac- 
teristic related to the coherence length of the bath as 
compared to the extension of the system, the equal cou- 
pling of each system oscillator to the latter might seem 
an arbitrary restriction. Imagine for example that each 
oscillator is at a different distance from the common heat 
bath, leading to an interaction 
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where the different oscillators feel a coupling of strength 
7i. The immediate consequence is that the effective cou- 
plings become 
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Though in this paper we will consider 7i = 1, the unbal- 
anced case would be solved following exactly the same 
procedure as we outline in the next section. 



III. NOISELESS SUBSYSTEMS AND 
ASYMPTOTIC PROPERTIES 

In this section we discuss the conditions to achieve 
noiseless subsystems where dissipation is avoided in one 
or two of the system's normal modes. The properties 
of our system are specified completely by the matrix 1-L 
appearing in Eq. ([2]): 
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Al2 ^2 A23 
Ai3 A23 CJ3 



(11) 



and we aim to derive the set of conditions for the system 
parameters leading to one or two normal modes decou- 
pled from the environment, i.e. whose effective coupling 
tin is zero. 

Let us consider a normal mode 5 with normal frequency 
Vts. The eigenvalue problem is expressed adequately by 
three equations: 
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components 



of the normal vector Ca = 



for 

(^15,^25,^35)^ with J^ij defined in Eqs.([5]). The con- 
dition for normal mode 5 to be non-dissipative (out of 
the bath' influence) leads to a constraint as follows 
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From Eqs. ([12]), (fT3|) and normalization condition we 
can obtain analytically C5, Vts with a further constraint 
for the system parameters. In other words, not all pa- 
rameter choices lead to NSs, but it is possible for some 
configurations of frequencies and couplings of the set of 
three oscillators (satisfying some constraint). 

The normal mode 5 in terms of the system parameters 
reads 



A13A12 + A23(^5 ~ ^1) 
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where c is the normalization constant. Applying Eq. (fT3|) 
we can obtain its eigenfrequency ^1 as: 
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where A = {uo\ 
titles S = {uj\ + 
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Eq. (fT6|) is one of our main results, and represents a hy- 
persurface in the d- dimensional parameters space (being 
d = {N + l)A^/2 = 6) whereby a normal mode is allowed 
to evolve freely and without dissipation. Such manifold is 
restricted to regions in which the normal mode frequency 
Qs is real and positive, and the normalization constant 
c is well defined. It's worth noting that looking only to 
non-dissipative modes (imposing the condition tvs = 0) 
has been crucial in order to solve analytically the above 
equations, otherwise we have to deal with complicated 
expressions involving third order equations (general nor- 
mal mode expression). 

When Eq. (p!6|) is fulfilled we obtain a NS composed 
by (at least) a single normal mode that is effectively un- 
coupled to the reservoir. This could be performed ar- 
tificially by tuning one of the d = 6 parameters of H, 
such as, for instance, the natural frequency of one oscil- 
lator. In experiments where local addressing is possible, 
such as ions confined to individual traps, this modifica- 
tion should be rather straightforward (see e.g. [45|). It 
should be stressed that noise models for ion traps typi- 
cally favor a SB interpretation in terms of fluctuating un- 
correlated surface dipoles |46j , though other microscopic 



models based on charge diffusion [47| in the electrode sur- 
face question whether the bath's correlation length could 
in fact be larger than the distance of the ion to the elec- 
trode. For the moment, this is an open problem. 

Configurations in which a NS consisting of two normal 
modes is produced can also be obtained analytically in 
the specific case of three oscillators. Indeed, when two 
normal modes uncouple from the environmental action, 
the third one must necessarily coincide with the center 
of mass of the system. Explicitly, the condition for the 
center of mass being a normal mode is: 
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with QcM = (1, 1, l)"^/\/3. The latter implication can 
be captured in the next two relations that have to be 
fulfilled simultaneously by the system parameters of the 
system: 



UJi = (jJ2^ A23 — Ai3 

CJ3 = a;2 + A12 — Ai3 



(18) 
(19) 



Since we want to remain in the domain of attractive po- 
tentials, we have to restrict ourselves to regions of the 
parameter space where ^'q^ = ^1 + ^3 — ^2 + ^^13 ^ 0- 

In order to see the power of the conditions (p!6|) and 
(|18|19p we will give in the following some examples of 
configurations in which a NS of one or two modes is pro- 
duced. We will consider for instance simpler situations 
in which the six-dimensional parameter space is reduced 
by assuming two of the three natural frequencies of oscil- 
lators to be equal, (cj = cji = 60^3 7^ CJ2) or in which two 
of the three couplings are equal (A = A12 = A23 7^ A13). 
This is sufficient to obtain some different scenarios ap- 
pearing in open an closed chains configurations, as is 
schematically shown in Figure [TJ 

Let's start from the case of two equal frequencies {uo = 
^1 = ^3 7^ ^2)- Then the quantities defined in Eq. (fT6|) 
are simply A = 0, S = cj^ — cj^ and 1Z = {— A13, A13 — 
A12 — A23}, the latter implying two different consistent 
solutions to Eq. ([16]): the first one {A12 = A23} and the 
second one {uj2 =002} with: 

~2 2 , 2Ai3(Ai2 + A23 - A13) — 2A12A23 .^^x 

^2 - ^ H T ^;^ tt: i^Uj 
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Both solutions can be simultaneously fulfilled too. In this 
case we have that: 



FIG. 1: (Color online) Different configurations for a chain 
of three coupled oscillators in which a NS of one (a-d) or 
two (e-f) normal modes is predicted. The tilde on parameters 
indicates a fixed value depending on the other non-tilded ones 
as described in the text. 



On the other hand, by assuming two equal couplings 
we have three different solutions: A = {0,A±} , where 
the first one is trivial and accounts for a separated pair of 
coupled oscillators and a non-coupled third. The second 
solution allows for situations Fig. [lb and Fig. [T]i, defined 
by: 



A± = Ai3 ± \ {UJ2 - ^l)(^2 - ^3) 



(22) 
= A~o} 



Finally for {uji = uJs A A13 =0 A X = uj"^ — uj^ 
we have a two- mode NS solution (Fig. [1^). 

The presence of one or two non dissipating normal 
modes prevents full thermalization of the system and 
leads to an asymptotic state whose features are analyzed 
in the following, focusing on entanglement and on quan- 
tum synchronization between the oscillators. 



A. Asymptotic Entanglement 

When a NS is enabled, decoherence can be prevented 
in the system leading to asymptotic entanglement that 
would be absent in the thermal state. As a measure of en- 
tanglement, we will use the well known logarithmic neg- 
ativity which is computable for bipartite Gaussian states 
as is our case: 
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E^ = max{0, — logi^_} 



(23) 



which satisfies Eqs. (jl8|19p and is thus a condition for a 
two- modes NS. These three situations correspond respec- 
tively to configurations in Fig. [1^ (A12 = A23), Fig. [TJd 
{ujI = cJ2^) and Fig. [if (A = Aq). It is worth noting 
that configuration in Fig. [1^ is valid also for the closed 
(Ai3 ^ 0) case, as well as Fig. [TJd for the open chain 
(when Ai3 = 0). 



where z^_ is the minimum symplectic eigenvalue of the 
partial transposed covariance matrix Vab^ correspond- 
ing to time reflection of one party. With the help of gen- 
eral expressions we can calculate analytically the asymp- 
totic entanglement when a NS is produced. Here we 
present our results for the open chain with equal cou- 
plings (A12 = A23 = A) and same frequencies for the 



external oscillators {uji = ujs = cj), when only one of the 
three normal modes is not subjected to dissipation (Fig. 
[1^) and when only one of them is dissipating (Fig. [1^) 
by imposing A = Aq. The details of the calculations are 
reported in Appendix [Al As initial condition for the nat- 
ural oscillators we choose a squeezed separable vacuum 
state given by: 
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where any other first order or second order moments are 
zero. 



1. One mode NS 

As a paradigmatic example of the case in which there is 
one frozen normal mode, let us consider the configuration 
given in Fig. [1^. As for the initial condition, uji = ujs = ^ 
in Eq. ([24|) and we will assume the same squeezing factor 
for the external pair, i.e. n = rs = r while the squeezing 
in the central oscillator r2 will be irrelevant. 

Normal modes coupled to the environment will reach 
an asymptotic thermal state whose variances are given 
in Eq. ()A3p while the uncoupled one evolves freely. The 
asymptotic covariance matrix of the external oscillators is 
obtained by expressing the second order moments of nat- 
ural oscillators in terms of the normal modes. Then we 
substitute respectively the asymptotic expressions corre- 
sponding to the frozen mode (not coupled to the bath) or 
the thermalized ones. This yields the following analytical 
expression for the entanglement: 

E^r = max{0, Eq + A£;(l + cos{2ujt))} (25) 

that is defined by a minimum value Eq and an oscillatory 
term with amplitude AE and frequency 2uj: 
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where re = {vq +rQ )/4 and the critical values are defined 
by the following expressions: 
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Coefficients ap and aq depend both on the bath's tem- 
perature and on the system parameters via the shapes 
and frequencies of the dissipative normal modes as can 
be seen in their definition in Eq. (jASp . Note that while 
decoupling of normal modes from the bath is tempera- 
ture independent, the amount of entanglement generated 
depends on it via the thermalized degrees of freedom. 

The presence of asymptotic entanglement between the 
external pair of oscillators in a symmetric chain (indepen- 
dent of the frequency of the central one, but depending 



on the temperature and initial squeezing) are illustrated 
in Figure [2l The minimum entanglement Eq is plotted 
both for low (left panel) and high temperatures (right 
panel) in the relevant squeezing ranges. Different regions 
are distinguished in the map and are labeled following 
Paz and Roncaglia notation in Ref. [23: sudden death 
is reached (SD), the asymptotic state consisting of an in- 
finite sequence of sudden death and revivals (SDR) and 
finally when non-zero entanglement is present at all times 
(no sudden death, NSD). 

An asymptotic entangled state with strictly Ej\f > 0, 
can be generated both when r > Tq" (> 2rc) or equiv- 
alently when r < Tq (< 2rc) with different origins. In 
the first case (r > Tq) the entanglement oscillates be- 
tween r — Tq" and r -{-r^ and the initial squeezing in the 
natural oscillators is employed as a resource to gener- 
ate an entangled state, while the bath contribution r^ 
is a source for its degradation. It's interesting to see 
that Tq" is strongly dependent on the bath's tempera- 
ture while the system parameters play a secondary role, 
only important at low temperatures. Indeed when tem- 
perature increases (T ^ uj) sudden death of entangle- 
ment can be only avoided by increasing r to be greater 

i log (4A2Tu; (^ + ^)) = i log(T) + ct. 

On the other hand the amplitude of the oscillations 
in this case is AE = 2rc that has a very weak de- 
pendence on temperature quickly reaching a constant 

value AE ^ I log (ou'^ (§- + ^) (^ + d^A when 
increasing the temperature. 

The second case {r < r^) is only present at low tem- 
peratures (of order O.luj). Here entanglement oscillates 
around r^ -\- r with amplitude 2r. This means that in- 
troducing no squeezing in the initial state leads to a con- 
stant entanglement at r^, while adding squeezing (a re- 
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FIG. 2: (Color online) Minimum entanglement Eq generated 
in the asymptotic limit between external oscillators of the 
chain in configuration in Fig. [T^ for low temperatures (left 
panel) and high temperatures (right panel). The different 
phases (SD, SDR and NSD) are bounded by the two critical 
values Tq (separating NSD phase from SDR phase) and — r^ 
(separating SDR from SD) that are represented by continuous 
white lines. Dotted line corresponds to 2rc and the dashed 
colored one with (rj — r^)/2. We have set 002 = 1.2(jJ and 
A = 0.6cj. 



source in the former case) makes entanglement tend to 
a SDR phase by widening its oscihatory amphtude. The 
fact that thermahzation can lead to entanglement at low 
temperatures is well known [50]. 

Finally, we can relate critical values Tq with the un- 
certainty induced by the environment in the virtual os- 
cillator q = {qi -\- (73) /a/2 position and momenta: 
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this reveals that when r^ > a squeezing in momentum 
is generated {Ap < ^) yielding entanglement as we have 
commented above. However note that we have never a 
minimum uncertainty state with r^ > r^ for all tempera- 
tures and physical regimes of system parameters. Indeed 
the uncertainty relation can be expressed for the virtual 
oscillator q as AqAp = ^e^o -^0 > i. The quantity re = 
(rg" +r^)/4 can be also related with virtual oscillator un- 
certainties in position and momenta as e~^'' = -^j^ < 1. 
In the left panel of Figure [2] we can see the two re- 
gions in which £^0 > (NSD phases): the big one at 
the left top corner corresponding to entanglement gen- 
eration by using the initial squeezing in the external os- 
cillators as a resource (once r > r^) and the small left 
bottom island that represents the environment yielding 
entanglement via the squeezing generated in the q coor- 
dinate when r < Tq . The SDR phase is centered around 
2rc (white dotted line) for low temperatures and their 
amplitude is given by the separation of the dashed col- 
ored line {r^ — r^)/2 from the zero squeezing axis. For 
temperatures greater than that for which r^ = (cross 
point between the dotted and dashed lines) they inter- 
change their roles acting now {r^ — r^)/2 (dashed col- 
ored Hue) as the center of the SDR region and 2rc as 
the amplitude. The SD phase is bounded by the quan- 



tity r^ - 4rc 



corresponding to the case in which 



^0 + 2AE < thus no entanglement is present in the 
asymptotic limit. For high temperatures (right panel) 
we can see how 2rc reaches a constant value while the 
pronounced curvature in the SDR region reveals that we 
can always obtain robust entanglement by increasing the 
squeezing parameter r logarithmically with temperature. 
We have to point out that our results resemble those 
obtained in Ref. [23 for two resonant harmonic oscil- 
lators. There a similar entanglement phases diagram 
has been found and the same two different mechanism 
for entanglement generation appear. In that context, 
both oscillations and the appearance of the low temper- 
atures NSD phase were attributed to non-markovian ef- 
fects, while here follow by considering a final asymptotic 
(Gibbsian) state for the normal modes coupled to the 
bath (that can be reproduced by a Markovian Lindb- 
land dynamics as is pointed in appendix [B]) . Presence of 
a third oscillator in the system allows for manipulation 
of the width of entanglement phases at low temperatures 
(specially the low squeezings NSD one) by tuning the free 



system parameters UJ2 and A. 



2. Two modes NS 

Let us now consider the case in which two modes be- 
come decoupled from the bath and in particular focusing 
on the configuration in Fig. [1^. This is indeed a symmet- 
ric open chain configuration as before, but now we have 
a special value of the couplings A = Aq leading to a larger 
NS. The calculation is similar to the previous one while 
now we have that only one of the normal modes thermal- 
izes and the other two have a free evolution decoupled 
from the bath. This leads to a less compact expression for 
asymptotic entanglement and more details are reported 
in Appendix [Al Still, a similar phase diagram for entan- 
glement can be found in this case by numerical evaluation 
of logarithmic negativity from Eq. (|A9p given in the ap- 
pendix [A) Results are shown in Fig. [3] in the same range 
of squeezing and temperatures as in the previous (one 
mode NS) case. For low temperatures (left panel) the 
low temperatures low squeezing NSD island of Fig. [2] 
that corresponds to the environment acting as a resource 
for entanglement generation disappears, since the bigger 
one expands to low squeezing. Degradation of resources 
by the environmental action here is not sufficient to pre- 
vent entanglement production even in the non-squeezed 
(r = 0) case for T < Tc since actually the mode e is also 
contributing to entanglement generation. On the other 
hand the entanglement phases shows the same structure 
for high temperatures (right panel) where the only differ- 
ence resides in the attenuated growth for entanglement 
when r increases (see color bars). 

Notice that all the expressions have been calculated in 
the limit of small gamma assuming a final Gibbs state 
for decohered eigenmodes and free evolution for nonde- 
cohered ones. Of course this situation can only be per- 
turbative, since for stronger coupling to the bath the 
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FIG. 3: (Color online) Minimum entanglement in the asymp- 
totic limit between external oscillators of the chain in the two 
modes NS configuration in Fig. [T^ for low temperatures (left 
panel) and high temperatures (right panel). The different 
phases (SD, SDR and NSD) are bounded by the continuous 
white lines obtained by numerical evaluation. We have set 
CJ2 = 1.2cj. 



eigenmodes become increasingly coupled among them, 
through second order processes mediated by the bath. 
This necessarily leads to decoherence of all eigenmodes, 
at a low rate though. 



B. Quantum synchronization 

In this section we analyze the dynamics of the system 
showing the existence of a parameter manifold where the 
oscillators oscillate in phase, synchronously, in spite of 
having different natural frequencies. Full dynamics for 
Gaussian states is characterized by mean values of po- 
sition, momenta and variances (second order moments). 
The possibility to have synchronization in this system 
is important for two reasons: (i) this phenomenon has 
been largely studied in non-linear systems but we show 
that for dissipation in a common bath, it can arise even 
among harmonic oscillators; (ii) few attempts have been 
done to extend it to the quantum regime and we show 
here that this phenomenon witnesses the presence of ro- 
bust quantum correlations, being actually a synchronous 
steady state accompanied by asymptotic entanglement. 

Spontaneous synchronization -where no external forc- 
ing is considered- has been recently studied in some quan- 
tum systems like optomechanical cells or nanomechanical 
resonators [33|, [3J, [5l|, [52| . Actually the phenomenon of 
mutual synchronization has been extended to the quan- 
tum regime considering a different scenario in Ref. [33| 
where two coupled harmonic oscillators with different fre- 
quencies and initially out of equilibrium are studied dur- 
ing their relaxation towards a thermal state. Synchro- 
nization was reported in first and second order moments 
during a long transient and accompanied by the robust 
preservation of quantum correlations (quantum discord 
[53|, [5J) between oscillators. Two oscillators dissipating 
in a common bath are actually preserving asymptotically 
their entanglement and retaining a larger ener gy than in 
the thermal state only if they are identical [29| . In this 
(symmetric) case they also evolve towards a synchronous 
asymptotic state [33|. When three elements are consid- 
ered, we have shown above that the symmetric chain can 
reach an asymptotic regime with entanglement between 
the external oscillators, independently on the frequency 
of the central one. Then asymptotic synchronization be- 
tween the external pair is also expected. Beyond this 
symmetric case, more interesting is the possibility offered 
by a chain to freeze all the oscillators out of the thermal 
state when their frequencies are all different, as discussed 
below. 

The long time dynamics of our system can be straight- 
forwardly calculated by assuming that normal modes 
which are coupled to the bath get thermalized, while un- 
coupled ones have a free evolution. This is sufficient to 
analyze the presence of synchronization in the asymp- 
totic state. In general quantum mutual synchronization 
appears always in one-mode NSs among natural oscilla- 
tors hnked by the non-dissipative mode, as long as they 



have an asymptotic dynamics with only one oscillatory 
contribution. Phase or antiphase synchronization at the 
non-dissipating normal frequency is possible in first or- 
der moments depending on the sign of their T matrix 
coefficients, while only in-phase synchronization occurs 
for second order moments at twice the frequency of first 
order ones. Let us illustrate it in some situations and 
compare with the time evolution of (gf) \li (Fig. [4|) 
when considering a simple Markovian dynamics in the 
weak coupling limit (see Appendix [B]) . 

Consider firstly the specific case of an open chain with 
equal couplings and frequencies in the external oscilla- 
tors (corresponding to situation in Fig. [1^). The form of 
the non-dissipative normal mode is C^^ = (1,0, — 1)^/\/2 
so synchronization will emerge only between external os- 
cillators in antiphase for the position and momenta at 
frequency Vl^ = ijj (the normal mode frequency) and for 
the second order moments (necessarily in-phase and at 
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FIG. 4: (Color online) Evolution of position variances for 
each oscillator in the open chain (see legend) for (a) con- 
figuration in Fig. [1^ where one mode NS is generated 
(cj = 1.3(jJ2 , A = 0A(jj2) synchronizing the external oscil- 
lators at 2(jO, (b) configuration in Fig. [1]: where a different 
mode NS is generated (cji = 1.2c<;2 , cjs = 1.Suj2 , A = OAUJ2) 
producing synchronization in all pairs of oscillators at 2Qe 
and (c) configuration in Fig. [1^ when a two modes NS is 
generated (cj = 1.3uj2 , A = Ao) and synchronization is lost. 
Bath parameters for the simulation are in all cases T — 10uj2, 
7 = 0.07c<;2 and A = 50cj2. 



2a;). The central oscillator instead decays into the ther- 
mal state, its initial oscillations being suppressed in the 
long time dynamics. This case is shown in Fig. [4^, where 
synchronization appears after a transient only for the ex- 
ternal oscillators of the open chain while the central os- 
cillator looses oscillation amplitude. 

In the latter case synchronization appears between 
identical unlinked (A 13) oscillators in a symmetric chain 
(Fig. [1^). More peculiar is the case in which all os- 
cillators have different frequencies and eventually cou- 
plings. In the case of Fig. [it we actually have that 
the non-dissipative mode involves all the three oscilla- 



tors Ce 



Ce(^: 



UJ2 , ^3 — ^2 ~ ^± 5 ^±) with 
Qe = v/<^2+^±- This can actually give rise to syn- 
chronous dynamics of all the oscillators, in spite of the 
difference in their natural frequencies. Since one of the 
components has different sign than the other two in Ce 
two of the oscillators first moments will synchronize in- 
phase between them and in antiphase with the third one. 
In Fig. IIJd a total synchronization is produced involving 
all three (different) oscillators, consistently with the fact 
that the non dissipative normal mode e involves all three 
oscillators. 

A different situation is produced when we have a two- 
modes NSs since two oscillatory contributions are present 
in the asymptotic limit of the natural oscillators. Here 
synchronization is only possible when the two normal 
modes frequencies are the same. An example is the open 
chain with a two-modes NS (Fig. [1^) where apart from 
the previous non-dissipative mode C^^, actually the mode 
Ce = (1,-2,1)^/^6 with frequency ft, = ^J^2Jl^^ 
does not dissipate either. In this case synchronization is 
destroyed by the presence of mode e and it can be only 
recovered when Vt, equals 1^^, i.e. in the trivial case of 
independent (A = 0) identical oscillators (u^ = cj). Lack 
of synchronization as well as a multimode oscillation are 
shown in Fig. [4h. 

The initial state employed for simulations is a squeezed 
separable vacuum state, where the squeezing parameters 
have been chosen to be different (t\ = 2,r2 = 2.5 and 
Tz = 3). In general, we have tried to avoid special ini- 
tial conditions that could have filtered just one normal 
mode into the dynamics. What we discussed is therefore 
the emergence of synchronization as a dynamical process 
when considering more general initial states, leading to 
robust conclusions. 

The scenarios here discussed allow to establish the ef- 
fect of having one or two modes NS in the configura- 
tions of open chains (Fig. [1^ and e). The same analysis 
can be extended to other cases where a different normal 
mode is uncoupled from the environment. For instance, 
the configurations in Fig. [lb and d admit only one non- 
dissipative normal mode that involves the three oscilla- 
tors, producing then a collective synchronization of the 
chain. 



IV. THERMALIZATION AND ROBUSTNESS 
OF QUANTUM CORRELATIONS 

Creation of NSs is a powerful tool to avoid decoherence 
and produce synchronized dynamics as we have seen in 
previous sections. However, the conditions leading to 
NSs are satisfied only in some parameter manifolds and 
it is relevant to see which is the effect of deviations from 
these couplings and detunings, that could also arise from 
the difficulty of experimental tuning. In this case, dissi- 
pation is present in all normal modes and the effective 
couplings of Eq. ([7j) are all different from zero. Thus 
a thermal state is finally reached in all the degrees of 
freedom. 

In absence of NS entanglement is lost after a finite 
time. Although the asymptotic state is simply the ther- 
mal (Gibbs) state, damping dynamics of the normal 
modes with different decoherence and relaxation time 
scales is present, producing a rich behaviour in which syn- 
chronization or high quantum correlations may emerge 
during a large transient before the final thermalization 
of the system. These effects have been reported recently 
in the case of two harmonic oscillators, where disparate 
decay rates between the two normal modes is produced 
for small deviations from the resonat case [33| . 

A dynamical description of the system in terms of a 
Markovian Master Equation (Appendix [Bj) reveals the 
central influence of the effective couplings in the relax- 
ation time scales of the different normal modes. In this 
context, the ratio between the two smallest decay rates, 
defined as: 






(32) 



provides important information about the dynamics of 
the system, being in fact one of the central figures in 
order to predict the robustness of correlations between 
oscillators or the emergence of synchronization as we will 
see in the following. In presence of disparate decay rates 
[R « 1), a large time interval appears between ther- 
malization of the two modes with largest damping co- 
efficients {strongly- damped modes) and thermalization of 
the mode with the smallest one {weakly- damped mode) . 
This produces the emergence, after a transient, of a long 
interval in which the weakly- damped mode is effectively 
the only one present in the dynamics, hence producing 
the synchronization between pairs of oscillators linked by 
this normal mode and the slow decay of quantum discord 
between these pairs. On the other hand, when the de- 
cay rates are similar {R ^ 1) the different modes are 
present for all times inhibiting synchronization, and the 
survival of correlations associated to one of the modes 
for long times is lost. These phenomena will be exempli- 
fied in the scenario of an open chain with equal couplings 
(A12 = A23 = A). 

In Figure [5] we represent R showing broad regions in 
which a weakly- damped mode exists (white regions) near 
the NSs manifolds corresponding to the configurations in 
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FIG. 5: (Color online) Map of R for A 
of the open chain frequencies. Dashed and dashed-dotted 
lines represent the non dissipative parameters manifolds of 
Figure [TJi and c respectively. The right panel is a zoom of the 
vicinities of the two-modes NS cross point in the left panel . 
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Fig. [T^ (dashed Hne),[Tt (dashed-dotted hyperbola) and 
[it (the crossing point). Out of these regions there is no 
separation of scales for the decay rates (blue regions) and 
all rates become progressively similar. We point out that 
the blue region wrapping the diagonal in Fig. [5] acts as 
boundary for the two white ones, since a different mode 
(with radically different shape) is weakly- damped in each 
white region. The coupling A is related to the position of 
the dashed-dotted hyperbola by Eq. ([22]) and the width 
of white regions, making them broader as A increases and 
tighter when it decreases. 



A. Quantum correlations 

Even if out of the NS conditions entanglement suffers 
a sudden death, following Ref. [3^ other quantumness 
indicators, such as quantum discord can remain robust in 
regions where disparate decay rates {R « 1) exist. By 
using an adapted measure of discord for Gaussian bipar- 
tite states [55|, [56| we observe the existence of a plateau in 
the dynamical evolution of discord between pairs of oscil- 
lators which are linked by a weakly- damped normal mode. 
More precisely, in the white region of Figure [3 close to 
the dashed-dotted hyperbola, the weakly- damped mode 
links the three natural oscillators, producing a plateau 
in the evolution of discord for all pairs. Moving to the 
tighter white region, close to the dashed diagonal line, 
the weakly- damped mode only involves the external oscil- 
lators pair of the open chain, leading to a slowly decaying 
discord only for this pair of oscillators. In blue regions 
no plateau is observed for discord, reaching for each pair 
the value corresponding to the thermal state in shorter 
times. 

Figure [6] shows time evolution of discord in logarith- 
mic scale for the three pairs of oscillators (see colors) for 
a selection of parameters close and far away from the 
dashed-dotted hyperbola (Fig. [6^ and [6)3 respectively). 
The initial condition has been taken to be a squeezed 



FIG. 6: (Color online) Time evolution of discord for the three 
pairs of oscillators in the open chain (see legend) for two dif- 
ferent regions of Figure [5l setting A = 0.4cji, cjs = I.6UJ2 
and changing cji. (a) Near the dashed-dotted hyperbola 
(cji = l.l(jj2) and (b) far away from it (cji = 1.9(jJ2). Bath 
parameters have been taken for the simulations T — 10uj2, 
7 = 0.07c(;2 and cutoff frequency A = 50uj2. 



separable vacuum state with same squeezing parameters 
as in Fig. [4] and will be kept for further simulations. A 
Gaussian filter has been employed to eliminate fast oscil- 
lations (original quantities are plotted in gray), in order 
to make it easier to identify the plateau characterizing 
discord robustness. 

As already seen for asymptotic entanglement, the ef- 
fect of increasing the bath's temperature is, in general, a 
degradation of quantum effects. It is therefore important 
to see how robustness of discord is a feature present also 
in hotter environment. The main effect when increasing 
T is that the final thermal state (by virtue of Eq. IB4p 
displays lower correlations so that the amount of discord 
that can be maintained in a robust way diminishes. In 
Figure [7| we show the evolution of discord for a pair of 
linked oscillators (1,2) for increasing T by factors 3 and 
6 (other parameters as in Fig. [6^): while the plateau is 
present for all temperatures and their (negative) slope is 
very similar, a lower amount of discord can be generated 
in the initial transient producing a shift of the curves to 
lower values (oscillations are increased by the logarith- 
mic scale of the plot). This degradation by temperature 
effects can be avoided by increasing the squeezing in the 
initial separable vacuum state like in the case of entan- 
glement presented in section [III Ai 



B. Synchronous thermalization 

With respect to the emergence of synchronization for 
pairs of oscillators when the NS is lost, we have to point 
out that when thermalizing the system reaches a station- 



10 




10- 



50 100 150 200 250 300 

UJot 



FIG. 7: (Color online) Time evolution of discord for a pair 
of linked oscillators (1,2) of the open chain for different bath 
temperatures (see legend). We have set A = 0.4u;|, UJ2, — 1.6c<;2 
and uj\ = l.lcc;2- The rest of bath parameters has been kept 
7 = 0.07c(;2 and A = 50uj2. 



ary state where oscillations are suppressed. We therefore 
restrict our analysis to a transient (which becomes longer 
the more we approach one of the NS conditions) where 
oscillations in the first and second order moments are still 
present in the dynamics. In this situation synchroniza- 
tion of first and second order moments can be estimated 
quantitatively by introducing a typical indicator, that is 
defined for two generic time-dependent functions h{t) and 
g{t) as: 



ChAt^^t) 



ir^'dfih- 



h){g-g) 



^/AhAg 



(33) 



ft+At 



where Ah = /;+''' dt'{h - hf and /i = ^ //^"^ dt'h{t'). 
When evolutions are phase or antiphase synchronized we 
will obtain \C\ ~ 1, while for very different dynamics we 
will obtain a value of C near to zero. 

Figure [8] shows the synchronization indicator C'^g2^^^2^ 

with position variances of (a) the external pair of oscil- 
lators i, j = 1,3 and (b) for i, j = 1,2 of the open chain 
with identical couplings and varying the external oscilla- 
tors frequencies (in the same range as in Fig. [5] ). We 
see immediately the high resemblance with the R map of 
Fig. [5] and some interesting differences induced by the 
shape of the normal modes. Effectively the external pair 
of oscillators (1,3) synchronizes (C ^ 1) in all regions 
where disparate decay rates (i? <C 1) are predicted since 
these two oscillators are linked by the weakly- damped nor- 
mal mode in these regions (Fig. [8^). As for the internal 
pair (1,2), it does depend on the weakly-damped mode 
in the vicinities of the dashed-dotted hyperbola where 
synchronization is actually found. On the other hand, 
near the diagonal the weakly- damped mode approximates 
to Cs = 1/a/2(1, 0, — 1) excluding the central oscillator 
from the induced collective motion. Consistently, the 1, 2 
pair (Fig. [8)3) does not synchronize for uji ~ uj2 (near the 
diagonal). 

We finally point out that, as expected, the synchroniza- 




1.0 1.2 1.4 1.6 1.8 2.0 1.0 1.2 1.4 1.6 1.8 2.0 

UJ1/UJ2 UJ1/UJ2 

FIG. 8: (Color online) Absolute value of the synchronization 
indicator |C(t,At)| for position variances ((g?)) for (a) the 
external pair of oscillators 1, 3 and (b) a linked pair 1,2. The 
synchronization factor is plotted at time t = min{tMAx, F^^} 
where Imax — 5000c(;2^^ (the maximum time used in the sim- 
ulations) in order to obtain a map in which oscillations were 



not yet suppressed. We have used At ■ 
bath parameters as in former figures. 



IScJo and the same 



tion frequency is that of the weakly- dam^ped mode (l^o) 
for the first order moments (position and momenta) and 
2Q.Q for the second order momenta. 



V. CONCLUSIONS 

Decoherence in an open quantum system can be 
avoided or reduced by tuning the system parameters in 
a common environment context. The shape of the in- 
teraction Hamiltonian between system and bath can be 
used in order to engineer the protection of some degrees 
of freedom from the environmental action. This analysis 
has been solved in the present work for a system of three 
coupled harmonic oscillators in contact with a bosonic 
bath in thermal equilibrium, developing the necessary 
general relations so as to obtain a NS composed by one 
or two non-dissipative normal modes. 

Different open and close chain configurations have been 
explored, highlighting the richer variety of NS configura- 
tions available when the dissipative system is extended 
from two to three harmonic oscillators. For a symmetric 
open chain (equal frequencies of the external pair of oscil- 
lators and same coupling) a closed analytical expression 
for the asymptotic entanglement (logarithmic negativity) 
has been derived observing the appearance of three dif- 
ferent phases depending on temperature and squeezing of 
the initial state (sudden death of entanglement, a infinite 
series of sudden death events and revivals and asymptotic 
robust entanglement). Sudden death of entanglement 
can be avoided for arbitrarily high bath temperatures 
by increasing the squeezing in the initial state for both 
one or two modes NS. Remarkably this critical squeezing 
in order to avoid sudden death depends logarithmically 
on temperature. 

Other dynamical effects such as the emergence of syn- 
chronization of mean values and variances have been ana- 
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lyzed in different situations by simply assuming the relax- 
ation to a thermal state of the normal modes coupled to 
the environment. Coherent oscillations appear when only 
a surviving normal mode is present in the dynamics, in- 
ducing synchronization in the natural oscillators that de- 
pend on it. Interestingly the parameter manifold leading 
to NSs include several not symmetric configurations: for 
instance an hyperbolic relation among frequencies can be 
satisfied for identical couplings in an open chain; in this 
case both asymptotic entanglement and synchronization 
are predicted even if all the oscillators natural frequen- 
cies are different (a possibility offered by a chain of three 
oscillators and absent in the case of two). A more general 
scenario where oscillators networks are considered with 
dissipation acting globally and locally has been discussed 
inRef.^. 

Furthermore an analysis of situations in which the NS 
conditions are not accomplished at all has been per- 
formed. Indeed, important properties can be present 
although when deviating from NS conditions entangle- 
ment does not survive: robust conservation of discord 
during a long transient dynamics and the emergence of 
synchronous oscillations are found before thermalization. 
These effects are interpreted in relation to disparate de- 
cay rates for the normal modes. As long as there is a 
weakly- damped mode surviving among several strongly- 
damped modes, effects such as robust discord and syn- 
chronization arise among the oscillators following this 
normal mode. 

While we have focused along this paper in specific 
cases of all the possible three oscillator configurations 
(in which calculations are greatly simplified) the strat- 
egy provided here is rather general and applies straight- 
forwardly to other choices of system parameters that pro- 
duce the decoupling of one or more normal modes. In fact 
this method can extend this engineering of the normal 
modes to more complicated systems such as disordered 
harmonic lattices or complex networks. 
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Appendix A: Analytical derivation of asymptotic 
entanglement 

All the information about bipartite quantum corre- 
lations for a Gaussian continuous-variable state is con- 
densed in its covariance matrix defined through the ten 
second order moments of q{A,B) and P{a,b) {^^ o^^ case 
first order moments are initially zero). This bipartite co- 
variance matrix defined for a system of two oscillators A 



and B, can be written as: 



Vab 



a 7 
7* 13 



(Al) 



where Q^,/3,7 are (2 x 2) blocks: a{P) contains the sec- 
ond order moments of oscillator subsystem A (B), and 
7 contains correlations of both subsystems. The mini- 
mum symplectic eigenvalue (of the covariance matrix cor- 
responding to the partially transposed density matrix), 
necessary to calculate the logarithmic negativity, is given 
by: 



yi(a ^b-2g- ./{^^^h^2gf 



As) (A2) 



with a = 4det(a), h = 4det(/3) , g = 4det(7) and s = 
16 det V . Normal modes coupled to the environment will 
reach in the asymptotic limit a thermal state, given by 
the Gibbs distribution e t /Tr(e t ). For a normal 
mode {k) we get: 



{Ql)th = 7;?rcoth( -^ 



2n 



2T 



(P,),, = -coth(- 



(A3) 



being Qk the corresponding normal mode frequency and 
T the reservoir temperature. On the other hand the un- 
coupled modes evolve freely. This means that the asymp- 
totic covariance matrix can be calculated by expressing 
all second order moments of natural oscillators in terms of 
the normal modes and then substituting the asymptotic 
expressions corresponding to free modes or thermalized 
ones. 

The covariance matrix in the asymptotic limit can be 
separated into three parts corresponding to the contri- 
butions of each normal mode. In terms of the blocks 
we can write: a = Yn^i^^Ai^i^ P = SLi-^-li^* ^^^ 
7 = 7^ = X]i=i J^AiJ^BiVi. For a non-dissipative normal 
mode (n) we have: 



K„ 



-dis 



(Ql) 



{{Qr.,Pr.}} 



{PD 



and for a dissipative one (fc) we get: 



Vdi. 



iQDth 







(Pi) 



th 



(A4) 



(A5) 



While elements in Vdiss are given by the expressions (|A3p 
those of Vno-diss are the ones corresponding to a free evo- 
lution of an harmonic oscillator. This analysis gives all 
the necessary elements in order to calculate asymptotic 
entanglement for pairs of oscillators in every particular 
situation in which one or two of the normal modes are 
un-coupled from the environmental action. 

Consider the specific case of an open chain (A13 = 0) in 
which we have two equal frequencies {ui = cJs = cj) and 



12 



two equal couplings (A12 = A23 = A 7^ 0) (Fig. [T^). In 
this case we get only one normal mode decoupled from 
the bath. In order to calculate the expression of the 
minimum symplectic eigenvalue we have to first calculate 
the elements of the three normal modes, that are shown 
here as vector columns: 




A 



; c± 



c± 



^l 



■ uo 



A 



where we have labeled the non-dissipative mode as 5 
and the other two modes as {±}. The correspond- 
ing normal modes frequencies are Q.s = ^ and Q.± = 
V(u;|+a;2)/2±v^, defining e = (^^^^)2 + 2X^, and 
c± is nothing but a normalization constant. We can now 
obtain all the terms appearing in V±. 

The initial condition given in Eq. ([24|) can be now re- 
written in terms of the non-dissipative normal mode as 



(Qm) 



-2r 



2uj 



, {pim 



u; 2r 
2^ 



and({Q5,P4(0)) 



0, and then their free evolution is given by: 



(Ql) 
(p!) 

{{Qs,Ps}) 



1 



2(jO 

|(e-2'-sin2(wi) + e2'-cos2(u;t)) 
2 sinh(2r) cos(wi) sin(wt) 



(A6) 



where we have already used that Qs = uj. By substituting 
in Vno-diss (Eq. IA4p we can now obtain the expressions 
of the determinants and introducing into equation (|A2p 
yields the following expression for the minimum symplec- 
tic eigenvalue: 



-{ty 



2A2 



Qo + Qi cos{2ujt) 



{Go + Gi cos{2ujt)y - 4apaQ (A7) 



which is an oscillatory function with frequency 2uj defined 
by the following functions: 



Go = {ctq - 

Gl = {CTQ 



- CTp) cosh(2r) 

- (jp) sinh(2r) 



and the dependence on bath temperature and shape of 
dissipative normal modes is given by: 



(jp 



2uj ^+ ^^^^ 1^ 2T 



^^c?_ coth 



2^. 



-ClCOth { 2T 



.2 



2^. 



-c_ coth 



\2t) 

(^) (A8) 



From Eq. (|A7p . we can obtain the minimum entangle- 
ment (given when t = (2n + l)^; n = 1,2,3, ...) and the 
maximum one (given when t = {n -\- 1 )^; n = 1,2,3,..) 
in order to recover Eq. ([25]) in section (jIII Ap of the main 
text with the proper definitions specified there. 

On the other hand, if we move to situation represented 
in Figure [1^ by fixing A = Aq (see Eq. [22]) we have that 



the normal modes transform into: 



C, = ^ 




a = ^ 



^CM 




being the last one, Ccm^ the only dissipative mode. The 
correspond ing norma l frequencies a re respect ively: fts = 



CJ, r^e = \/2ijJ2 



P' and Vt 



CM 



^2^ 



Naturally 



we have to restrict ourselves to the regime 2CJ3 > a; > 
uj\l2 in order for these quantities to be real and positive. 
Keeping the same initial condition as in the previous 
case, we have that nothing changes in the expression of 
the free evolution of mode (5 of Eq. ([A6p while the free 
evolution of mode e is given by: 



{pD 

({Qe,P.}) 



_ 2a£2+Ci£^2r 

2uJ^UJ'2^ 



f±^e2-sin^(l]et) + f^e-2-cos2(l]et) 



6... e-2-sin2(l]et) + ^^e2-cos2(l]et) 

g2r _ (2a.+a;.)n. ^-2r^ ^0^(^1,1)^1X1(^1,1) 



f 2ou2-\-uJ 2r 



We have assumed the same squeezing parameter r in the 
central oscillator of the chain (note that in the previous 
case the initial state of the central oscillator is not rele- 
vant and we have not specified it). Following the same 
procedure as above, we calculate an expression for the 
minimum symplectic eigenvalue. It is worth noting that 
in this case we have two contributions to the determi- 
nants of the free type Vno-diss (Eq. IA4p corresponding 
to the two non dissipative modes, and only a dissipative 
one Vdiss (Eq. IA5p corresponding to the center of mass 
mode. 

The minimum symplectic eigenvalue yields: 



2zy_(t)2 



Ao^Aiit) 



(A9) 



-^{Ao^Ai{t))^-Bo-Bi{t) 



where we have defined the following quantities in order 
to simplify the expression. The constant terms: 



^0 

Bo 



cosh(2r) (4(aQ + dp) + J7+(1]^ + uj^)) 



GAapCTQ 



81uJijU2 

and the oscillating terms: 



V ^e 



Ai{t) = 4cos(2a;t)sinh(2r)(crQ - ap) + 
+ J+ cos{2ujt) sinh(2r)(^^ + cj^) + 
+ JL cos{2net) cosh(2r)(l]^ - uj^) - 

i^e^Oj)^ 



{Q,-cu)^' 



- j-cos{2{ne-oo)t)smh{2r)- 

- J- cos{2{ne ^ uj)t) smh{2r) 

where different frequencies coming from the two non- 
dissipative modes are present. We have used J± = 



13 



and the two new bath de- 



12UJ \ rig UJUJ2 

pendent functions are now simply given by the contribu- 
tion of the center of mass mode's thermal state: 



a 



CFp 



^Q 



CM 

Quo 



coth 



2T 



-cothf^^ 



(AlO) 



of frequencies characterized by an spectral density J{f2). 
For simplicity it has been considered to be Ohmic with 
a sharp cutoff J(i7) = — i? 0(A — 1]), where Q{x) is the 
Heaviside step function, A is the largest frequency present 
in the environment (cutoff frequency) and 7 is a constant 
quantifying the strength of system-environment interac- 
tion (thus in the weak-coupling we have always that 7 
<C r^i Vz G {1,...,A^}). This assumptions leads to the 
following definitions of the Master Equation coefficients: 



Appendix B: Master Equation 

Assuming the general framework provided in section 
HH we can postulate a simple approach to the dynamical 
evolution in the weak-coupling limit between system and 
environment. By using the general Born and Markov ap- 
proximations as well as initial product state, a Markovian 
Master Equation (MME) for the reduced density matrix 
of the open system can be obtained [4^ in the normal 
modes basis. The resulting equation is not of the Lindb- 
land form, thus complete positivity (CP) is not guaran- 
teed [53. This issue can be solved in two different ways, 
either considering a rotating- wave-approximation (RWA) 
in the interaction Hamiltonian (j4]) XiXj -^ aiO/- ^ h.c. or 
performing a strong-type RWA in the non-Lindbladian 
master equation by eliminating oscillatory terms of the 
form e±*(^i±^i)^ that appear in the interaction picture. 
The latter is the one we will pursue. The advantages of 
this method not only reside in obtaining a Master Equa- 
tion in Lindbland form (thus CP), but also in that dy- 
namical evolution can be solved analytically. However 
an exhaustive analysis in the case of two harmonic os- 
cillators shows a very well agreement between results us- 
ing the original non-Lindbladian master equation and the 
strong RWA here used [33j . 

The MME for the evolution of the reduced density ma- 
trix for a common bath in the strong RWA is then: 



-i[Hs,p{t)]- 

2 E r« ([Q«. {^". PW}] - [Pn, {Qn, Pit)}]) 



dp{t) 
dt 



+ Dn [[Qn,[Qn,Pm " ^ [^n, [P„, p(t)]] ) (Bl) 

Here F^ and Dn are constant coefficients (by virtue of the 
Markov approximation) accounting for the damping and 
diffusion effects respectively. Note that under this ap- 
proximation, each normal mode is dissipating separately 
to the bath, i.e. they have independent decay rates. The 
bath has been considered to be in a thermal (Gibbs) state 
at temperature T, and to be composed by a continuum 



L r) — 1^, 



Dr, 



2 TT J{Qn) 



7^n 



nr 



(B2) 



^Ji^n) coth (^1=7 f^l^n coth ( -^ 



2T 



2T 



where we also assume 1^^ < A Vz G {1, ..., N}. 

As we are interested in classical and quantum correla- 
tions of the system oscillators, a description for the evo- 
lution of the first and second order moments is necessary. 
This can be written in a simple form as: 

R = MR + N. 

where R is a column vector containing the M = {2N -h 
1)N independent second order moments of the normal 
modes. The matrix A^ condenses all the information 
about their dynamical evolution and N determines the 
stationary values for long times (when R = 0). The 
dynamics of R can be solved in terms of the eigenvalues 
ofA^: 



{^ij} = {- 



±i\Qi±Qj\} • i<j 



(B3) 



where the i = j eigenvalues determine the evolution of 
(Qf), (P^) and {{Qi.Pi}), while the ones with i ^ j de- 
termine that of (QiQj), {PiPj) ^^d {{Qi^Pj}). Note that 
by virtue of eqs. (|B2|) and (|B3P the decay of the normal 
modes is entirely governed by the effective couplings men- 
tioned above, thus differences in their magnitude produce 
disparate temporal scales for the dissipation and diffusion 
of normal modes. 

The stationary state is found to be (Roo = A^~"'"N): 



(Q')c 



D, 



2F,a2 



1 , f^i 

coth — 

2^, V 2T 



D 



^i 



(i^^)^ = ^ = :iicoth^ 



2F 



^i 



2T 



being all the other second order moments equal to zero. 
Note that these expressions for the asymptotic limit re- 
cover the thermal state of the system at the bath tem- 
perature T given by the Gibbs distribution in Eqs. (|A3p . 
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